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Radiation pressure and absorption in AGN: results from a complete 
unbiased sample from Swift 
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ABSTRACT 

Outward radiation pressure can exceed the inward gravitational pull on gas clouds in the 
neighbourhood of a luminous Active Galactic Nucleus (AGN). This creates a forbidden region 
for long-lived dusty clouds in the observed columnn density - Eddington fraction plane. (The 
Eddington fraction AEdd is the ratio of the bolometric luminosity of an AGN to the Eddington 
limit for its black hole mass.) The Swift/BAT catalogue is the most complete hard X-ray 
selected sample of AGN and has 97 low redshift AGN with measured column densities TYh 
and inferred black hole masses. Eddington fractions for the sources have been obtained using 
recent bolometric corrections and the sources have been plotted on the N-r — AEdd plane. 
Only one source lies in the forbidden region and it has a large value of Ah due to an ionized 
warm absorber, for which radiation pressure is reduced. The effective Eddington limit for the 
source population indicates that the high column density clouds in the more luminous objects 
lie within the inner few pc, where the central black hole provides at least half the mass. Our 
result shows that radiation pressure does affect the presence of gas clouds in the inner galaxy 
bulge. We discuss briefly how the Nn — AEdd plane may evolve to higher redshift, when 
feedback due to radiation pressure may have been strong. 
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^ 1 1 INTRODUCTION 

I_h ' Evidence for a relatively tight relation between the mass of the 
central black hole of a galaxy and the mass or velocity disper- 



sion of the surrounding stellar bulg e (Kormendy & Richstone 1995 



Mago rrian etafll l998. Geb hardt et al. 2000, iFerrarese et alJuOOl] 



sure from a central Active Galactic Nucleus (AGN) on surrounding 
dusty gas in a hard X-ray, flux-limited sample of AGN. 



2 THE EFFECTIVE EDDINGTON LIMIT 



Tremaine et al. 2002) points to a coupling or feedback of the black 
hole on the galaxy. The accretion power produced by the growth 
of the black hole exceeds the binding energy of the galaxy by 
one to two orders of magnitude so some form of feedback is ex- 
pected. What is un c lear however is how it works. Is it the energy 
dSilk & Reeslll998l, lHaehnelt et aTl 1 1998h or momentum (Fabian 
ll999l : lFabian et aljl2002klKindl2003l ; lMurrav et alj|2005h supplied 
by the central engine that are most effective? Winds, radiation pres- 
sure and jets have all been invoked and it is likely that they all play 
a role, although perhaps at different phases and in different objects. 

Direct evidence of feedback is elusive and difficult to obtain. A 
quasar at the centre of a galaxy will have little effect on the stars of 
the surrounding bulge and only act on the gas. Catching and recog- 
nising an object where the gas is being expelled may not be easy. 
Here we provide clear evidence for the effect of radiation pres- 



iLaor & Drainel (1 1993b . IScoville & Normanl dl995h and 
Murra vet alj J2005h have shown that the effective Eddington 
limit, when the outward radiation pressure on gas exceeds the 
inward gravitational pull, can be much lower for dusty gas than 
for ionized dust-free gas, as assumed for the standard derivation. 
A reduction factor of 500 is obtained for an AGN spectral energy 
distribution (SEP), peaking in t he UV and assuming a Galactic 
dust-to-gas ratio jFabian et al.1 [2008). X-ray ionization keeps 
the drift velocity of the dust grains relative to the gas low. The 
UV radiation is rapidly used up by a relatively small column 
density of gas and dust (equivalent hydrogen column density of 
Nh ~ 10 21 cm~ 2 ) and the effective Eddington limit drops off 
approximately linearly, reaching unity at column densities where it 
goes Compton thick around A^h ~ 10 24 cm -2 . This means that 
an AGN which has a luminosity one to two orders of magnitude 
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below that standard Eddington limit may nevertheless exceed the 
effective Eddington limit for substantial columns of dusty gas. 

The actual interaction of the accelerated face with the rest of 
an absorbing cloud will be complicated and likely Rayleigh-Taylor 
unstable, with the details depending on magnetic fields, geome- 
try and on how the force is transmitted to the whole cloud. The 
net result however is that long-lived stable clouds are not expected 
to survive in a regime where the effective Eddington limit is ex- 
ceeded. This manifests as a forbidden reg ion in the column den- 
sity - Eddington fraction plane for AGN (Fabian et ah 2008|). By 
Eddington fraction, AEdd , we mean the ratio of the bolometric lu- 
minosity of the source to the standard Eddington limit. An AGN 
with AEdd = 1/500 is at the Effective Eddington limit for Galac- 
tic abundance dusty gas with Ah ~ 10 21 cm -2 . This is shown in 
Fig. 1. This limit has been calculated assuming an average SED typ- 
ic al for low Eddington ratio A GN (the thick black curve in Fig. 13 
of I Vasudevan & Fabianl 120071) : at high Eddington ratios the sub- 
stantially different SED shape may marginally reduce the effective 
Eddington limit . 

The above argument applies only where the black hole dom- 
inates the mass locally. For outer gas hundreds of pc to kpc from 
the black hole, the AEdd boundary increases in proportion to the 
enclosed mass of stars. Outer (few kpc) dust lanes in a galaxy 
bulge or disk are unlikely to have very high column densities (the 
total gas mass at radius rkpc is otherwise very high, M ga s ~ 
10 10 /iV23J- 2 M Q , where the column density is 10 23 A 23 cm" 2 and 
/ is the covering fraction), so we suggest a lower boundary to the 



forbidden region of Ah ~ 
upper limit to N H of 10 21 - 10 2 
gala xies by [ Winter et al. (2008a). 



. This agrees with the 
found for Seyfert 1-1.2 



iFabian et al . (2008) plotted several samples on the Ah — Asdd 
plane. These were a composite low-redshift sample of 23 objects 
from Swift and BeppoSAX, a sample of 77 sources from the Chan- 
dra Deep Field South and 13 objects from the Lockman Hole. The 
results provided tentative evidence for sources avoiding the for- 
bidden region. Column densities are derived from X-ray spectra 
and Eddingto n fractions from X-ray lum inosity (and a bolomet- 
ric correction. Ivasudevan & Fabian 2007) and importantly an esti- 
mate of the black hole mass. This was mainly done by combining 
K-band ma gnitudes of the AGN host with the Mbh — K corre- 
lation of iMarconi & Hunt! <2O03h . IWinter et all <2008ah find that, 
before application of a bolometric correction, the distributions of 
ix/^Edd for absorbed and unabsorbed sources are significantly 
different. Here we include the bolometric correction to recover the 
full Eddington ratio. We present our analysis of 97 sources from 
the 9- month Swift BAT catalogue of AGN, w ith BH mass esti- 
mates dTueller et ai]|2008L IWinter et ai]|2008al) . which is the most 
complete and uniform, hard X-ray flux-limited, sample of AGN. 
These sources lie above the Galactic Plane (|6| > 15 deg) and have 
2MASS data for mass determination. A hard X-ray sample is im- 
portant since it will be relatively insensitive to absorption (provided 
that the sources are not too Compton thick). 



3 RESULTS FROM THE SWIFT/BAT 9-MONTH AGN 
CATALOGUE 

In Fig. [2] we present the column density Ah against Eddington 
ratio for the 9-month Swif t BAT catalogue of AGN. The column 
densities are presented by Wint er et alj d2008al) and were obtained 
from fits to X-ray spectra from a variety of X-ray missions in- 
cluding ASCA, Suzaku, Chandra, XMM-Newton and the X-ray tele- 
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Figure 1. Nn — AEdd plane showing the forbidden region where dusty 
clouds close to the black hole see the AGN as being effectively above the 
Eddington limit (blue solid line). Outer dust lanes can obscure over much 
of the plane at low column densities (black horizontal line) and long lived 
absorbing clouds can only occur in the upper left region, where the AGN is 
below the effective Eddington limit. The effective Eddington limit depends 
on the dust grain abundance: the solid blue line shows the effective Ed- 
dington limit for a standard ISM grain abundance, and the dot-dashed and 
dotted blue lines are for grain abundances of 0.3 and 0.1 of ISM abundance 
respectively. 



scope (X RT) on Swift. The XRT fits were performed by Wint er et al.l 
( 2008a) themselves, whereas the other fits were taken from the liter- 
aturo The Eddington ratios are calculated u sing the standard for - 
mulaA Ed d = Li,„i/1.38 x 10 38 (Mbh /M g ) Winter etahl d2008al) 
estimate the black hole mass Mbh using the A/bh — ibuigc rela- 
tion, using the K-band magnitudes from the 2MASS All-Sky Sur- 
vey to provide bu lge luminosities. The A/bh — -^buige relation from 
INovaketalJ d2006» is employed by these authors to calculate their 
BH masses, which are presented in Table 4 of their work. We then 
scale up the 2-10 keV luminosity to a bolometric luminosity us- 
ing a hard X-ra y bolometric correction d ependent on the X-ray 
spectral shape. Ivasudevan & Fabianl d2007l) find that the bolometric 
correction (and optical/UV/X-ray SE D shape in gene r al) de pends 
on Eddington ratio, and the study of She mmer et alj d2006l) (and 
IShemmer et alj |2008 subsequently) finds a correlation between the 
X-ray photon index, r, and Eddington rati o. This latter correlation 
is not reproduced bv lWinter et alj d2008al) , possibly due to the in- 
trinsic scatter in individual AGN accretion rates and photon indices, 
or a systematic deviation of lower luminosity AGN from the corre- 
lation. However, they do find F correla tes with AEdd in ind ividual 
objects varying between observations d Winter et al. 2008b), lend- 
ing credence to the possibility that spectral shape is linked to ac- 
cretion rate, albeit in a non-straightforward fashion. Based on these 
findings, we adopt a bolometric correction of 19 for hard spectrum 



1 The source represented in Fig. 2 by an empty circle near (0.1, 22.8 
shown after updating its parameters using archival XMM-Newton data. 
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Figure 2. Absorbing column density against Eddington ratio for the Si- 
month Swift/BAT catalogue of AGN. Circl es represent objects for which 
X-ray fits were taken from the literature by Winter et al. 1 2008a), whereas 
squares represent those for which Swift XRT spectra were used. Filled 
circles and squares are those reported to have complex X-ray absorption 
and empty circles an d squares are those re ported to have simple absorbed 
power-law spectra, in [Winter et alj J2008ah . The objects with arrows point- 
ing upward s in the lower region of the plot did not have any column density 
reported in Winter et al. 1 2008j3) and we show them with a nominal upper 
limiting column density of 10 20 cm -2 . The dashed line shows a factor of 
2 increase in the effective Eddington limit due to the mass of intervening 
stars. All other conventions are as given in Fig.fT] 



objects (r < 1.9) and 55 for soft spectrum objects (r > 1.9) and 
thus determine the Eddington ratios. 

We note that the AGN from the 9-month BAT catalogue pri- 
marily lie in the region expected for long-lived absorption or dust 
lanes. There are a few objects which lie above the effective Edding- 
ton limit for dusty gas, but we note that this limit has hitherto been 
calculated assuming that the only inward gravitational force of im- 
portance is that from the central black hole. In reality, there may be 
a significant gravitational force from the stars located inward from 
the dusty gas clouds of interest. This would increase the effective 
Eddington limit by an amount proportional to the total enclosed 
mass, since the increased gravitational force is now able to balance 
a proportionately higher outward force due to radiation. 

We estimate the importance of this effect as follows. If we 
enclose a mass equal to twice the central BH mass, the effective 
Eddington limit is doubled. Such a limit would allow the major- 
ity of the objects to now lie within the effective Eddington limit, 
with the exception of one (MCG-05-23-16). We can also comment 
on the constraints that such a modified effective Eddington limit 
w ould place on t he loc ation of the dusty gas clouds. The studies 
of ISchodel etafl J2007l) reveal that, in the case of our Galaxy, the 
enclosed mass doubles at a scale of a few parsecs from the nu- 
cl eus. If nuclear s tar clusters often accompany AGN as the study 
of lSefhet alj d2008h seems to suggest, then an effective Eddington 
limit that accounts for twice the BH mass would still comfortably 
locate the dusty gas clouds close to the nucleus. 



4 DUST-TO-GAS RATIO AND THE EFFECTIVE 
EDDINGTON LIMIT 

We note that the dust-to-gas ratio in AG N may deviate signifi - 
cantly from the Galactic dust-to-gas ratio dMaiolino et alj 1200 ll). 
The c alculations of the effective Eddington limit in iFabian et al.1 
(2008) assume the Galactic value, via the assumption of inter- 
stellar ^jriedhrm_^ISM)gmn abundance in their CLOUDY simula- 
tions. |MaioMno|eral] ([2001) suggest that the ratio of optical extinc- 
tion to X-ray column density A*h is systematically lower in AGN, 
and here we account for this effect by running further CLOUDY 
simulations with dust grain abundances reduced by factors 0.3 and 
0.1 of the standard ISM abundance (Fig.[TJ. The effect of this alter- 
ation is signficant at low A^h, where the effective Eddington limit is 
inversely propotional to the dust to gas ratio. At high column densi- 
ties, the effective Eddington limit with the reduced dust-to-gas ra- 
tio converges with that for Galactic dust abundances, since atomic 
absorption in the gas dominates. However, since the change in ef- 
fective Eddington limit at low column densities coincides with the 
regime where dust lanes could be responsible for the measured A*h , 
this modification has little bearing on the shape of the 'forbidden' 
region in A*h — A E dd space. We note that the effective Eddington 
limit will also depend on the mass of the black hole, since the peak 
of the thermal black body in the UV depends on BH mass, but do 
not address this further here. The average SED used in the CLOUDY 
simulati ons is calculated from eigh t AGN SEDs with a range of BH 
masses dVasudevan &Fa bian 2007) so is probably a representative 
choice of SED. 



5 DISCUSSION 

Sources in the most complete sample of hard X-ray selected AGN 
available avoid the forbidden region in the A/h — A E dd plane. This is 
unlikely to be due to an observational selection effect since sources 
there should be brighter. The single source which lies in the forbid- 
den region is MCG-05 -23-016, which h as a complex spectrum with 
warm absorbers (e.g. lReeves et al.l2007l) and so not necessarily ab- 
sorbed by dusty gas clouds. Our results provide strong evidence 
that radiation pressure from the central AGN is the main agent. 

The good agreement between the distribution of sources and 
our allowed zones suggests that our assumptions are resonable. The 
radiation pressure seems to be transmitted to the whole cloud, or 
cloud complex, perhaps by magnetic fields. The Eddington frac- 
tions must be accurate to a factor of a few dWinter et al. 2008a 
suggest a factor of 3 uncertainty) and therefore the bolometric lu- 
minosities also. Our result (Fig. 2) provides a simple explanation 
for some of the otherwise complicated behaviour that has been 
report ed in the relatio nship between absorption and luminosity in 
AGN (Hasinger 2008 and references therein). 

There can be matter out of our line of sight in the low column 
density, high Eddington fraction sources but that matter will need 
to avoid the forbidden zone. For example it can reside in a Compton 
thick torus, as envisaged in many unification models for AGN. 

We note that the soft X-ray bright Narrow-Line Seyfert 1 
galaxies, which are generally considered to have high Eddington 
fractions do not show up in a hard X-ray selecte d sample ow- 
ing to their steep X-ray spectra Mali zia et alj d2008l) . They do not 
have strong absorption so would appear at the lower right of the 
A^h — A E dd plane. 

We do not know the evolution of absorbing clouds in these 
galaxies, nor the evolution of the central black holes. Thus we can- 
not say whether the lowly-absorbed AGN have prevented clouds 
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from accumulating or have blown away an existing population of 
clouds. The UV radiation from the AGN cannot vary greatly on the 
local dynamical timescale, 10 4 r pc (T 1 ^ ) fl yr where r pc and aioo are 
the radius in pc and velocity dispersion in units of 100 kms , in 
order to obtain the agreement we see. Objects can of course appear 
in the forbidden zone for timescales much shorter than this. The 
absence of such objects argues that the luminosity of these AGN is 
either fairly steady or the duty cycle of high luminosity episodes is 
small. 

What should in future become possible with large samples 
from XMM, Chandra and IXO is to study the evolution of the 
population of AGN. The distribution of AGN in the TVh — Asdd 
plane seen as a function of redshift should reveal in which direc- 
tion they evolve. The mass doubling timescale for the black hole is 
~4x 10 7 Ag^ d yr, assuming an accretion efficiency of 0.1. Since 
many of the objects in Fig. 2 have AEdd < 0.01, then they are 
growing on a timescale > 4 Gyr, which is slow. How cold gas ac- 
cumulates within the host galaxy bulge and settles near its centre 
will determine the fuelling of the AGN and central star formation 
rates. Cold gas typically resid es for about 10 yr before forming 
stars in low AEdd galaxies (e.g. lKennicutt|[l998l) . so star formation 
occurs faster than the black hole mass evolves. The rates at which 
the mass of the black hole and host galaxy bulge are changing are 
however both low. 

The evolution of quasars means that high Eddington rate ob- 
jects with short doubling times are common around redshifts of 
two and we can expect much higher densities of objects near the 
effective Eddington limit. The black hole mass is changing faster 
than gas is turned into stars and the effect of radiation pressure on 
the expulsion of gas can have dramatic effects on t he mass reser- 
voir for growing st ars as well as the black hole (see iFabianll 1 999L 
iMurrav et alj|2005l and lFabian et alj|2008l for possible scenarios). 
The locus of the effective Eddington limit in the Nh — AEdd plane 
specifies the radius of the absorbing gas within the host galaxy and 
thus the masses of gas involved. Although we cannot deduce the 
evolution of an individual object from its position in the plane, how 
the population evolves may be attainable. 

Our results show direct evidence for a coupling of the radia- 
tion from AGN and the surrounding gas. 
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